Asymptotic Bound-state Model for Feshbach Resonances 
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We present an Asymptotic Bound-state Model which can be used to accurately describe all Fesh- 
bach resonance positions and widths in a two-body system. With this model we determine the 
coupled bound states of a particular two-body system. The model is based on analytic properties of 
the two-body Hamiltonian, and on asymptotic properties of uncoupled bound states in the interac- 
tion potentials. In its most simple version, the only necessary parameters are the least bound state 
energies and actual potentials are not used. The complexity of the model can be stepwise increased 
by introducing threshold effects, multiple vibrational levels and additional potential parameters. 



The model is extensively tested on the Li- 
system are presented. 

PACS numbers: 



K system and additional calculations on the K- Rb 



I. INTRODUCTION 

The field of ultracold atomic gases has been rapidly 
growing during the past decades. One of the main sources 
of growth is the large degree of tunability to employ ul- 
tracold gases as model quantum systems [ij [2] . In par- 
ticular the strength of the two-body interaction param- 
eter, captured by the scattering-length a, can be tuned 
over many orders of magnitude. A quantum system can 
be made repulsive (a > 0), attractive (a < 0), non- 
interacting (a = 0) or strongly interacting (|a| ^ oo) in a 
continuous manner by means of Feshbach resonances [3 . 
These resonances are induced by external fields: magnet- 
ically induced Feshbach resonances are conveniently used 
for alkali-metal atoms, while optically induced Feshbach 
resonances seem more promising for e.g. alkaline-earth 
atoms. In this paper we consider magnetically induced 
resonances only. 

Feshbach resonances depend crucially on the existence 
of an internal atomic structure, which can be modified 
by external fields. For alkali-metal atoms, this structure 
is initiated by the hyperfine interaction, which can be 
energetically modified by a magnetic field via the Zee- 
man interaction. For a given initial spin state, its col- 
lision threshold and its two-body bound states depend 
in general differently on the magnetic field. A Feshbach 
resonance occurs when the threshold becomes degenerate 
with a bound state. Accurate knowledge of the Feshbach 
resonance structure is crucial for experiments. 

The two-body system has to be solved to obtain the 
bound state solutions. Since the interactions have both 
orbital and spin degrees of freedom, this results in a set of 
radially coupled Schrodinger equations in the spin basis. 
The set of equations is referred to as Coupled Channels 
equations [4 , and can be solved numerically. Quite of- 
ten it is far from trivial to obtain reliable predictions 
for the two-body problem, due to several reasons: the 
ab-initio interaction potentials are usually not accurate 
enough to describe ultracold collisions. Therefore these 
potentials have to be modelled by adding and modify- 



ing potential parameters. A full calculation for all spin 
combinations and all potential variations is very time- 
consuming. Moreover, one can easily overlook some fea- 
tures of the bound state spectrum due to numerical issues 
such as grid sizes and numerical accuracy. This is also 
due to a lack of insight of the general resonance structure, 
which is often not obvious from the numerical results. 

Given the above, there is certainly a need for fast and 
simple models to predict and describe Feshbach reso- 
nances, which allow for a detailed insight in the resonance 
structure. In the last decade various simple models have 
been developed for ultracold collisions [5-7 , which vary 
significantly in terms of complexity, accuracy and appli- 
cability. In all these models the radial equation plays a 
central role in describing the Feshbach resonances. 

In this Paper we present in detail the Asymptotic 
Bound-state Model (ABM). This model, briefly intro- 
duced in Ref. [8 , and extended in Ref. [9 was success- 
fully applied to the Fermi- Fermi mixture of ^Li and ^^K. 
In Ref. [8 the observed loss features were assigned to 
13 Feshbach resonances with high accuracy, and the ob- 
tained parameters served as an input to a full coupled 
channels analysis. The ABM builds on an earlier model 
by Moerdijk et al. |T0] for homonuclear systems, which 
was also applied by Stan et al. [11 for heteronuclear sys- 
tems. This earlier model neglects the mixing of singlet 
and triplet states, therefore allowing the use of uncou- 
pled orbital and spin states. In the ABM we make use 
of the radial singlet and triplet eigenstates and include 
the coupling between them. This crucial improvement 
makes the whole approach in principle exact, and it al- 
lows for a high degree of accuracy given a limited number 
of parameters. 

We show how we can systematically extend the most 
simple version of ABM to predict the width of the Fesh- 
bach resonances by including threshold behavior. Addi- 
tionally we allow for the inclusion of multiple vibrational 
levels and parameter for the spatial wavefunction overlap. 
The fact that ABM is computationally light provides the 
possibility to map out the available Feshbach resonance 



positions and widths for a certain system, as has been 
shown in Ref. [9]. Throughout the paper we wih use 
the ^Li-^^K mixture as a model system to ihustrate ah 
introduced aspects. Additionahy, we present ABM cal- 
culations on the ^^K-^^Rb mixture to demonstrate its 
validity on a more complex system, comparing it with 
accurate coupled channel calculations [12j. The case of 
metastable helium atoms where each atom has an elec- 
tron spin of s = 1 and the interaction occurs through 
singlet, triplet and quintet interaction potentials we dis- 
cuss elsewhere [13 . 

In the following we describe the ABM (Sec. |ll|) and 
various methods to obtain the required input parame- 



ters. In Sec. Ill the ABM is applied to the three physi- 



cal systems and in Sec. [IV|we introduce the coupling to 
the open channel to predict the width of Feshbach res- 
onances. In Section |V| we summarize our findings and 
comment on further extensions of the model. 



II. ASYMPTOTIC BOUND-STATE MODEL 

In this section we give a detailed description of the 
asymptotic bound state model. In Section p[I A| we start 
with a general overview of the model which is described 



in more detail in the subsequent sections II B to II F 



A. Overview 

In the ABM we consider two atoms, a and /3, in their 
electronic ground-state. To search for Feshbach reso- 
nances we use the effective Hamiltonian 14i 



n = n 



rel 



-w- 
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Here H^^^ = P^/2/i + V describes the relative motion of 
the atoms in the center of mass frame: the first term is 
the relative kinetic energy, with /i the reduced mass, the 
second term the effective interaction potential V. The 
Hamiltonian W^^ stands for the internal energy of the 
two atoms. 

We will represent W^^ by the hyperfine and Zeeman 
contributions to the internal energy (Section II B ). There- 
fore, W^^ is diagonal in the Breit-Rabi pair basis {|Qf/3)} 
with eigen-energies Eaf3 and typically dependent on the 
magnetic field B. The internal states |a/3) in combination 
with the quantum number / for the angular momentum of 
the relative motion define the scattering channels {a/3^ I). 

Because the effective potential V is in general not di- 
agonal in the pair basis {|Qf/3)}, the internal states of the 
atoms can change in collisions. To include the coupling 
of the channels by V, we transform from the pair basis 
to a spin basis {\cr)} in w hich Ti^^^ is diagonal. We will 
restrict ourselves (Section [ll C ) to effective potentials V 
which are diagonal in S, the quantum number of the total 
electron spin S = SQ,-hS/3, where s^ and s^ are the elec- 
tron spins of the colliding atoms. The effective potential 



can thus be written as V(r) = "^^ \S)Vs{r){S\, where 
r is the interatomic separation. The examples discussed 
in this paper are alkali atoms {s = 1/2) which lead to a 
decomposition in singlet {S = 0) and triplet potentials 
{S = l). 

The eigenstates of W^^ (bound-states and scattering 
states) are solutions of the Schrodinger equations for 
given value of /, using effective potentials Vq jr) in which 
the centrifugal forces are included (Section II C). Since 
the effective potentials are central interactions, a sep- 
aration of variables can be performed to describe the 
wavefunction as a product of a radial and angular part, 
1^) = |^)|y^^). The ABM solves the Schrodinger equa- 
tion for the Hamiltonian (IT]) starting from a restricted set 
of (typically just a few) discrete eigenstates IV^^Ol^m^) ^^ 
H^®\ using their binding energies e^^ as free parameters. 
The continuum states are not used in the model. The 
set {IV^^O} corresponds to the bound-state wavefunctions 
i^u^i^) = {A^v^) i^ ^^^ effective potentials V|(r), with 
V and / being the vibrational and rotational quantum 
numbers, respectively. 

The ABM solutions are obtained by diagonalization 
of the Hamiltonian (fTl) using the restricted set of bound 
states {\'4^v^)\(y)} ' This is equivalent to solving the secular 
equation 

det|(ef-E5)(^.Za,.'ZV' + (V^f |V^5')(^I'^^"V0I =0, (2) 

where we have used the ort honor mality of |F^^). The 
roots El) represent the eigenvalues of 1-L which are shifted 
with respect to the hare levels e^^ due to the presence of 
the coupling term H^"*. The roots £^5 will be accurate 
as long as the influence of the continuum solutions is 
small. Since the bound-state wavefunctions il)y\r) are 
orthonormal for a given value of S and /, the Franck- 
Condon factors are (V^^^lV^f/) = Sw' • The eigenstates of 
1-L define bound states in the system of coupled channels. 

We define the entrance channel {a/3, Oo ^^ ^^^ internal 
states I a/3) for which we want to find Feshbach resonances 
with a given angular momentum state of / = 0, 1, • • •. The 
energy E^JB) of two free atoms at rest in the entrance 
channel defines the threshold energy, which separates the 
continuum of scattering states from the discrete set of 
bound states. In the ABM we define 1-L relative to this 
energy. With this convention the threshold energy al- 
ways corresponds to £ = 0, irrespective of the magnetic 
field. Further, we consider only entrance channels that 
are stable against spin exchange relaxation. 

Since in the ABM we only consider bound states, and 
therefore are in the regime £^ < 0, all channels are ener- 
getically closed. Collisions in the entrance channel would 
have a collision energy of £ > and the entrance chan- 
nel would be energetically open, i.e. the atoms are not 
bound and have a finite probability of reaching r = 00 
in this channel. Although all channels are closed in the 
ABM we will refer to the entrance channel as the open 
channel, anticipating on the inclusion of threshold effects 
in Section ITVl 



In the Sections III BIIII Fl we discuss the ABM in its sim- 
plest form, where level broadening by coupling to the con- 
tinuum is neglected [8]. In this approximation, Feshbach 
resonances are predicted for magnetic fields Bq where a 
bound level crosses the threshold, £^5 = /J^rei{B — Bq) 
with /j^rei = dE})/dB\B=Bo^ and where coupling to the 
continuum is in principle allowed by conservation of the 
angular momentum. To determine the crossings the di- 
agonalization Q has to be carried out as a function of 
magnetic field. 

The procedure becomes particularly simple when the 
coupling strength W^^ is small compared to the separa- 
tion of the ro- vibrational levels in the various potentials 
V^{r) because in this case the basis set can be restricted 
to only the least bound level in each of the potentials 
y|(r). In this case the set of levels {e^^} reduces to a 
small number, {e*^^}, with |5q, — 5/3I < S < 5q, + 5/3. In the 
case of the alkalis only two levels, e^^ and e^^, are relevant 
for each value of / . Further, as will be shown in Section 
|IIF[ the least bound states have the long-range behavior 
of asymptotically bound states, which makes it possible to 
estimate the value of Franck-Condon factors ('^^^l'^^/^ 
without detailed knowledge of the short-range behavior 
of the potentials V|(r). This reduces the diagonalization 
([2]) to the diagonalization of a spin Hamiltonian. Treat- 
ing the {e^^} as fitting parameters, their values can be 
determined by comparison with a hand full of experimen- 
tally observed resonances. Once these {e^^} are known, 
the position of all Feshbach resonances associated with 
these levels can be predicted. As this procedure does not 
involve numerical solution of the Schrodinger equation 
for the relative motion it provides an enormous simplifi- 
cation over coupled channels calculations. 



In Section [IVl we turn to the extended version of the 
ABM in which also the coupling to the open channel 
is taken into account. The presence of such a coupling 
gives rise to a shift A of the uncoupled levels and above 
threshold to a broadening F [9 . The width of a Fesh- 
bach resonance is related to the lifetime r = h/T of the 
bound-state above threshold and provides a measure for 
the coupling to the continuum. In magnetic field units 
the width AB is related to the scattering length by the 
expression [10, 



B. Internal energy 

To describe the internal energy of the colliding atoms 
we restrict the atomic Hamiltonian to the hyperfine and 
Zeeman interactions 



w 






i • s + (7eS - 7ii) • B, 



(4) 
(5) 



where s and i are the electron and nuclear spins respec- 
tively, 7e and 7^ are their respective gyromagnetic ratios, 
tthf is the hyperfine energy and B is the externally ap- 
plied magnetic field. The hyperfine interaction couples 
the electron and nuclear spin which add to a total angular 
momentum f = s + i. In Fig.[l]the well known Breit-Rabi 
diagrams of ^Li and ^^K are shown, the curves correspond 
to the eigenvalues of 1-L^. The one- atom hyperfine states 
are labeled \fmf), although / is only a good quantum 
number in the absence of an external magnetic field. 
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Figure 1: The single atom hyperfine diagrams for Li and 
^°K. The curves correspond to the eigenvalues of H^ and are 
labeled by the zero field quantum numbers \fmf). 



By labeling the colliding atoms with a and /3, the two- 
body internal Hamiltonian becomes W^^ = H^ + H^ 
and the spin state of the colliding pair can be described 
in the Breit-Rabi pair basis \aP) = \fa'^fo,-, fjs'^fp) = 
I /, ^/ ) a ^ I /, ^/ ) (3 ' The corresponding energy of two free 
atoms at rest defines the 5-dependent threshold energy 
introduced in Section III A[ 



a{B) = abg [1 



AB 
B-Bo 



(3) 



where a^^ is the background scattering length. Interest- 
ingly, the width AB can also be determined with the 
same restricted basis set {\i^^^)}, which does not include 
continuum states. In Section IIVI this is shown for the 
simplest case where only a single level is resonant and 
the resonance width can be found from the coupling of 
two bound-state levels: the resonant level and the least- 
bound level in the entrance channel. 



C. Relative Hamiltonian 

The bound eigenstates of Ti^^^ play a central role in 
the determination of the coupled bound states H respon- 
sible for the Feshbach resonances. The relative Hamil- 
tonian includes the effective interaction V resulting from 
all Coulomb interactions between the nuclei and electrons 
of both atoms [33] . This effective interaction is isotropic 
and depends only on the quantum number S associated 
with the total electron spin. For these central poten- 
tials the two-body solutions will depend on the orbital 
quantum number /, but not on its projection mi. In the 



absence of any anisotropic interaction both / and rui are 
good quantum numbers of H^^^ and H. 

We specify the ABM basis states as {li^^^)^)}. Here 
the spin basis states \a) = \SMsfial-(^(3) are determined 
by the spin quantum number S and the magnetic quan- 
tum numbers M^, //q,, and /i/3 of the S, Iq, and i/3 op- 
erators, respectively. The sum Mp = Ms + Ma + M/3 is 
conserved by the Hamiltonian H. This hmits the number 
of spin states which have to be included in the set \a). 
The bound-state wavefunctions il^y\r) for the singlet and 
triplet potentials, characterized by the vibrational and 
rotational quantum numbers v and /, satisfy the reduced 
radial wave equation of H^®^ for specific values of 5* and 
I 



2/i dr'^ 



+^iw 



SLlSU 



r.\r) = et^r.'ir). 



(6) 



Here V|(r) = ]/s(r) + /(/ + l)fi^/(2/ir^) represents the in- 
teraction potentials including the centripetal forces. The 
corresponding binding energies are given by e^K In this 
paper we mainly focus on heteronuclear systems, how- 
ever, the ABM works equally well for homonuclear sys- 
tems. In the latter case one would rather use a sym- 
metrized spin basis \a) = \SMsIMi), where / is the to- 
tal nuclear spin and Mj is the magnetic quantum number 
for I = icK + i^ as described in Ref. ^lOj . 



D. Diagonalization of 1-L 

In the ABM basis {\i^^^)\o')} the Zeeman term Hz is 
diagonal with 

E^ = {(7\nzW) = KleMs - Tc/ic - 1^^^)B (7) 

the Zeeman energy of state |cr). As the orbital angular 
momentum is conserved, we can solve Eq. S separately 
for every / subspace. Since the set {\ipf,^)\cr)\\^ orthonor- 
mal the secular equation takes (for a given value of /) the 
form 

det|(ef + Ef -Eb),5,,,,V'+<'5{'^|^hfk')| =0, (8) 

where r]^'^, = {ip^ \ip^f ) are Franck- Condon factors be- 
tween the different S states, which are numbers in the 



s,s' 



s,s 



Re- 



range < \r]^y I < 1 for 5 7^ S'^ and rj^y 
peating the procedure as a function of magnetic field the 
energy level diagram of all bound states in the system of 
coupled channels is obtained. 

It is instructive to separate the hyperfine contribution 
into two parts, Hhf = Hj[f + Hj^f , where 



'^hf - 



^(s«±s;3)-i.±^(s.±s^)-i^. (9) 



Because H^^ conserves S^ it couples the ABM states only 
within the singlet and triplet manifolds. The term H^^ 
is off-diagonal in the ABM basis, hence couples singlet 



to triplet. Accordingly, also the hyperfine term in the 
secular equation separates into two parts 

(10) 



Note that the second term of Eq. (10) is zero unless 
S y>^ S'. This term was neglected in the models of 
Refs. [TOlITT]. This is a good approximation if no Fesh- 
bach resonances occur near magnetic fields where the en- 
ergy difference between singlet and triplet levels is on the 
order of the hyperfine energy. However, for a generic case 
this term cannot be neglected. 

To demonstrate the procedure of identification of Fesh- 
bach resonances we show in Fig. |2] the ABM solutions for 
a simple fictitious homonuclear system with S = 1 and 
/ = 2 for an entrance channel with Mp = Ms + Mj = 
and / = 0. The example has the spin structure of ^Li 
but we use ABM parameters, e^, e^ and 77^^, with values 
chosen for convenience of illustration. The field depen- 
dence of threshold energy of the entrance channel E^^ 
is shown here explicitly (dashed line). The energies E^ 
(solid lines) are labeled by their high field quantum num- 
bers \SMsIMi) and the binding energies in the singlet 
and triplet potentials are chosen to be e^ = —10 and 
e^ = —5. The avoided crossings around 5 = are caused 
by the hyperfine interaction and is proportional to ahf; 
the avoided crossing between the singlet and triplet lev- 
els is proportional to the wavefunction overlap rf^. Four 
s-wave Feshbach resonances occur indicated at the cross- 
ings 1, 2 and 3 (double resonance). The resonances at 1 
and 2 are mostly determined by the triplet binding en- 
ergy e^ and the resonances at 3 by the singlet binding 
energy e^. 



E. Free parameters 

The free parameters of the ABM are the binding en- 
ergies e^'^ and the Franck- Condon factors r]^'^, . These 
parameters can be obtained in a variety of manners. Here 
we discuss three methods, two of which will be demon- 
strated in Sect. Ill and the third in Ref. [T3l. 



First, if the scattering potentials V|(r) are very well 
known, the bound state wavefunctions of the vibrational 
levels can be obtained by solving equation (p| for e^^ < 
0. Numerical values of the Franck- Condon factors follow 
from the obtained eigenfunctions. This method is very 
accurate and can be extended to deeply bound levels, 
however accurate model potentials are only available for 
a limited number of systems. 

A second method can be used when the potentials are 
not very well or only partially known. For large inter- 
atomic distances the potentials can be parameterized by 
the dispersion potential 



V(r) 



a 



^6 • 



(11) 



Since this expression is not correct for short distances. 
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Figure 2: ABM calculation for a fictitious homonuclear sys- 
tem with 5* = 1 and 1 — 2 for an entrance channel with 
Mf = Ms + Mi = and / = 0. The threshold energy of the 
entrance channel E^f^ is shown here explicitly as the dashed 
line. The energies Eh (solid lines) are labeled by their high 
field quantum numbers \SMsIMi). The binding energies of 
the least bound states in the singlet and triplet potentials are 
chosen to be e° = — 10 and e^ = —5. The avoided crossing 
around B — is proportional to the hyperfine interaction auf 
and those between the singlet and triplet levels to the wave- 
function overlap 7701 and the hyperfine interaction auf- Four 
Feshbach resonances occur indicated at the crossings 1, 2 and 
3 (double resonance). 



we account for the inaccurate inner part of the poten- 
tial by a boundary condition based on the accumulated 
phase method [15 . This boundary condition has a one- 
to-one relationship to the interspecies s-wave singlet and 
triplet scattering lengths. This approach requires only 
three input parameters: the van der Waals Cq coefficient 
and the singlet {as) and triplet {ar) scattering lengths. 
For an accurate description involving deeper bound states 
the accumulated phase boundary condition can be made 
more accurate by including additional parameters [15 . 

The third method to obtain the free parameters is by 
direct comparison of ABM predictions with experimen- 
tally observed Feshbach resonances, for instance obtained 
in a search for loss features in an ultracold atomic gas. 
A loss feature spectrum can be obtained by measuring, 
as a function of magnetic field, the remaining number 
of atoms after a certain holding time at fixed magnetic 
field. The ABM parameters follow by a fitting procedure 
yielding the best match of the predicted threshold cross- 
ing fields with the observed loss feature spectrum. We 
applied this method in Ref. [8], where it has proven to 
be very powerful for rapid assignment of Feshbach reso- 
nances in the ^Li-^^K mixture due to the small compu- 
tational time required to diagonalize a spin hamiltonian. 

The number of fit parameters is determined by the 



number of bound states which have to be considered. 
Depending on the atomic species and the magnetic field, 
only a selected number of vibrational levels e^^ have to 
be taken into account. This number can be estimated by 
considering the maximum energy range involved. An up- 
per bound results from comparing the maximum dissoci- 
ation energy of the least bound vibrational level D"" with 
the maximum internal energy of the atom pair £^int,max- 
The maximum dissociation energy of the z^-th vibrational 
level can be estimated semi-classically [T6j 



^* ^ 



vCh 



v/^ 



1/2^ 



1/6 



(12) 



where C = 2 [r(l + l/6)/r(l/2 + 1/6)] ^ 3.434 where v 
is counted from the dissociation limit, i.e. z^ = 1 is the 
least bound state. The maximum internal energy is given 
by the sum of the hyperfine splitting of each of the two 
atoms at zero field, the maximum Zeeman energy for the 
free atom pair and the maximum Zeeman energy for the 
molecule 

^int,max ^ K/ + ^^hf + 2(5, + Sp)gs^BB, (13) 



where E\ 



hf 



= a 



hf 



(^Q;,/3 + <5a,/3) and we have neglected the 



nuclear Zeeman effect. Comparing equations (12) and 



( 13 ) gives us an expression for the number of vibrational 



levels Nj^ which have to be considered 
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1/2^1/6 
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1/3 n 
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(14) 



where \x\ denotes the smallest integer not less than the 
argument x. The maximum possible magnetic field Bmax 
required to encounter a Feshbach resonance can be esti- 
mated from Eq. [12] by neglecting the hyperfine energy 
as 



Br, 



D* 



{Sa + Sf3)gslJ^B ' 



(15) 



If the hyperfine energy is comparable or larger than the 
vibrational level splitting I^* the expression for B^ax 
overestimates the maximum field of the lowest field Fesh- 
bach resonance. 



F. Asymptotic bound states 

The most crucial ABM parameters are the binding en- 
ergies e^K However, for accurate predictions of the Fesh- 
bach resonance positions also the Franck- Condon factors 
have to be accurate. For weakly bound states these fac- 
tors are mainly determined by the difference in binding 
energy of the overlapping states, rather than by the po- 
tential shape. Therefore good approximations can be 
made with little knowledge of the scattering potential. 



For very- weakly bound states the outer classical turn- 
ing point Tc is found at distances Vc^ vq; i.e., far beyond 
the Van der Waals radius of the interaction potential 



To 



1 / 2/iC6 

2 V h^ 



1/4 



(16) 



These states are called halo states^7\. Because in this 
case most of the probability density of the wavefunc- 
tion is found outside the outer classical turning point, 
these states can be described by a zero-range potential 
with a wavefunction of the type ip{r) ~ e~'^^, where 
K, = (— 2/ie/fi^)^/^ is the wavevector corresponding to a 
bound state with binding energy e. The Franck-Condon 
factor of two halo states with wave vectors f<io and /^^i is 
given by 



{^'\^') 






(17) 



This approximation is valid for binding energies |e| <C 

The calculation of the Franck-Condon factors can be 
extended to deeper bound states by including the disper- 
sive van der Waals tail. For distances r :^ rx, where 
rx is the exchange radius, the potential is well described 
by Eq. (11) and the Franck-Condon factors can be cal- 



culated by numerically solving the Schrodinger equation 
(|6| for the Van der Waals potential (11) on the interval 
rx < r < oo [18 . The exchange radius rx is defined as 
the distance where the Van der Waals interaction equals 
the exchange interaction. This method can be used for 
asymptotic bound states, which we define by the condi- 
tion re > rX' If even deeper bound states, with Tc < rx, 
have to be taken into account, the potential can be ex- 
tended by including the exchange interaction [19 , or by 
using full model potentials. 

To illustrate the high degree of accuracy achieved by 
using asymptotic bound states we calculate the Franck- 
Condon factor for a contact potential (halo states), a 
van der Waals potential (asymptotic bound states) and a 
full model potential including short range behavior, de- 
rived from Refs. [20, 21 . Figure [s] shows the Franck- 
Condon factor r]il for ^Li-^^K calculated numerically 
for the model potential, van der Waals potential, and 
analytically using equation (17). The van der Waals 
coefficient used is Cq = 2322EhaQ 21], where Eh = 
4.35974 X 10-^^ J and ao = 0.05291772 nm. The value of 
r]il has been plotted as a function of the triplet binding 
energy e^ for three different values of the singlet binding 
energy e^. It is clear that the contact potential is only ap- 
plicable for e/h < lOOMHz, hence only for systems with 
resonant scattering in the singlet and triplet channels. 

The approximation based on the Cq potential yields 
good agreement down to binding energies of e/h < 
40GHz, which is much more than the maximum possi- 
ble vibrational level splitting of the least bound states 
{D'^/h = 8.2GHz), hence is sufficient to describe Fesh- 
bach resonances originating for the least bound vibra- 



tional level. The black circle indicates the actual Franck- 
Condon factor for the least bound state of ^Li-^^K. For 
the contact, van der Waals and model potentials we find 
r]^l = 0.991, r]'^l = 0.981 and r]^l = 0.979 respectively. 
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Figure 3: (Color online) The Franck-Condon factor r]2\/ for 
the least bound states of the ^Li-^'^K system, calculated as 
a function of the triplet binding energy e^ for three dif- 
ferent values of e^/h = 7.16 MHz, e^/h = 716 MHz and 
e^/h = 7.16 X 10^ MHz. r]^^^, is calculated for the model 
potential (dashed blue), the — Ce/r ^ potential (solid red) and 
the contact potential, equation (17) (dash-dotted green). The 
black circle indicates the actual value for the least bound state 
of ^Li-^°K {e^/h = 716 MHz and e^/h = 425 MHz). The 
nodes in ri^\/ correspond approximately to the appearance 
of deeper lying vibrational states, i.e. for e^ /h > 10^ MHz, 



III. APPLICATION TO VARIOUS SYSTEMS 

In this section we demonstrate the versatility of the 
ABM by applying it to two different systems using the 



different approaches as discussed in Section HE 



A. ^Li 



^K 



Both ^Li and ^^K have electron spin s — 1/2, there- 
fore the total electron spin can be singlet 5* = or triplet 
S* = 1. We intend to describe all loss features observed in 
Ref. [8 . Since all these features were observed for mag- 
netic fields below 300 G we find that, by use of Eq. (14), 



it is sufficient to take into account only the least bound 
state {u = 1) of the singlet and triplet potentials. This 
reduces the number of fit parameters to e^ 



SI _ 0,1 



lA 



and 
Subsequently, we calculate the rotational shifts by 
parameterizing the / > bound state energies with the 
aid of model potentials [34 as described by [20l |2T] . 
This allows us to reduce the number of binding energies 
to be considered to only two: e{ = e^ and e{ = e^. 
We now turn to the Franck-Condon factor rj^l of the two 



|1/2,+1/2>Li + |9/2,5/2>K |1/2,+1/2>Li + |9/2,-7/2>k 5=1 
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Figure 4: (Color online) The energies of all the coupled 
bound states for ^Li -^"^K with total spin Mp — ±3. The 
black solid line indicates the threshold energy of the en- 
trance channel |l/2,+l/2)Li |9/2,+5/2)k for 5 < and 
|l/2,+l/2)Li |9/2, -7/2)k for 5 > 0. The grey area rep- 
resents the scattering continuum and the (colored) lines in- 
dicate the coupled bound states. Feshbach resonances occur 
when a bound state crosses the threshold energy. The color 
scheme indicates the admixture of singlet and triplet contri- 
butions in the bound states obtained from the eigenstates of 
the Hamiltonian (IlJ. The strong admixture near the thres- 
hold crossings at 5 ^ 150 G demonstrate the importance 
of the singlet-triplet mixing in describing Feshbach resonance 
positions accurately. Since in these calculations the coupled 
bound states are not coupled to the open channel, they exist 
even for energies above the threshold. 



bound states. As discussed in Section IIIFI its value is 
^11 — 0-979 and can be taken along in the calculation 
or approximated as unity. We first consider the case of 
r]W = 1, this reduces the total number of fit parame- 
ters to only two. We fit the positions of the threshold 
crossings to the 13 observed loss features reported in 
Ref. [8 by minimizing the x^ value while varying e^ and 
e^. We obtain optimal values of e^/h = 716(15) MHz and 
€^/h = 425(5) MHz, where the error bars indicate one 
standard deviation. In Fig.[4| the threshold and spectrum 
of coupled bound states with Mp = +3 (—3) is shown 
for positive (negative) magnetic field values. The color 
scheme indicates the admixture of singlet and triplet con- 
tributions in the bound states. Feshbach resonances will 
occur at magnetic fields where the energy of the coupled 
bound states and the scattering threshold match. The 
strong admixture of singlet and triplet contributions at 
the threshold crossings emphasizes the importance of in- 
cluding the singlet-triplet mixing term Ti^^ in the Hamil- 
tonian. All 13 calculated resonance positions have good 
agreement with the coupled channel calculations as de- 
scribed in Ref. [8 , verifying that the ABM yields a good 
description of the threshold behavior of the ^Li— ^^K sys- 
tem for the studied field values. 

We repeat the x^ fitting procedure now including the 
numerical value of the overlap. The value of rj^l for both 



the 5-wave and p-wave bound states are calculated nu- 
merically while varying e^ and e^. This fit results in a 
slightly larger x^ value with corresponding increased dis- 
crepancies in the resonance positions. However, all of the 
calculated resonance positions are within the experimen- 
tal widths of the loss features. Therefore, the analysis 
with rjil = 1 and rj^l = 0.979 can be safely considered to 
yield the same results within the experimental accuracy. 



B. ^"K -^^Rb 

We now turn to the ^^K -^^Rb mixture to demonstrate 
the application of the ABM to a system including mul- 
tiple (three) vibrational levels in each potential and the 
corresponding non-trivial values for the Franck-Condon 
factors. We consider s-wave (/ = 0) resonances. Al- 
though accurate K-Rb scattering potentials are known 
[23 , we choose to use the accumulated phase method as 
discussed in Section \Ll E| using only three ABM parame- 
ters to demonstrate the accuracy of the ABM for a more 
complex system like ^^K -^^Rb. 

We solve the reduced radial wave equation |6l for 
Vs{r) = -Cq/t^ and the continuum state E = /i^/2/i 
in the limit A: ^ 0. We obtain the accumulated phase 
boundary condition at Tin = 18 a^ from the bound- 
ary condition at r ^ 00 using the asymptotic s-wave 
scattering phase shift Sq = arctan(— /ca), where a is the 
known singlet or triplet scattering length. Subsequently, 
we obtain binding energies for the three last bound states 
of the singlet and triplet potential by solving the same 
equation (|6| but now using the accumulated phase at 
r = rin and '^(r ^ 00) = as boundary conditions. We 
numerically calculate the Franck-Condon factors by nor- 
malizing the wavefunctions for r^n < r < roo^ thereby 
neglecting the wavefunction in the inner part of the po- 
tential (0 < r < Tin)- This approximation becomes less 
valid for more deeply bound states. We use as input 
parameters Cq = 421 A E^a^ [22], as = —111.5 a^ and 
ar = —215.6 ao [23 . Figure [s] shows the spectrum of 
bound states with respect to the threshold energy for 
the spin mixture of |/, tti/) = |9/2, — 9/2)k and |l,l)Rb 
states. The red curves indicate the ABM results and 
the blue curves correspond to full coupled channel calcu- 
lations [24 . The agreement between the two models is 
satisfactory, especially for the weakest bound states close 
to the threshold. A conceptually different analysis of the 
K-Rb system using also only three input parameters has 
been performed by Hanna, et al. [7|. 



IV. FESHBACH RESONANCE WIDTH 

A. Overview 

The asymptotic bound state model has been used so 
far to determine the position of the Feshbach resonances 
but not their width. As is well known from standard 
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Figure 5: (Color online) The bound state spectrum for ^"^K 
-^^Rb for Mf — —7/2 plotted with respect to the threshold 
energy ^K,Rb of the |9/2, -9/2)k + |l,l)Rb mixture. Solid 
red lines are ABM calculations and the blue dashed lines are 
numerical coupled channels calculations. Good agreement be- 
tween the two calculations is found in particular for the weak- 
est bound levels. 



Feshbach theory, the width of 5-wave resonances depends 
on the coupling between the resonant level and the con- 
tinuum [25l |26] . For resonances with / > the width is 
determined by a physically different process, namely tun- 
neling through the centrifugal barrier. Here we discuss 
only the width of s-wave resonances. We determine the 
resonance width by analyzing the shift of the resonant 
level close to threshold due to the coupling with the least 
bound state of the open channel. This is possible using 
again the restricted basis set of bound states {|'0^O|cr)}, 
introduced in Section|n[ The possibility to obtain the res- 
onance width by analyzing the shift is plausible because 
near a resonance the scattering behavior in the zero en- 
ergy limit is closely related to the threshold behavior of 
the bound-state. To reveal the coupling as contained in 
the ABM approximation we partition the total Hilbert 
space of the Hamiltonian ([l]) into two orthogonal sub- 
spaces V and Q. The states of the open channels are 
located in V space, those of the closed channels in Q 
space [To]. This splits the Hamiltonian 1-L in four parts 
(cf. Section[lVB|; H = Hpp+Hpq+Hqp+Hqq. Here 
T-Lpp and 1-Lqq describe the system within each subspace 
and 1-Lpq{= Hqp) describe the coupling between the V 
and Q spaces, thus providing a measure for the coupling 
between the open channels in V space and the closed 
channels in Q space. 

The scattering channels are defined by the Breit-Rabi 
pair states |a/3) = |/a^^/« •> f^'^fp)- I^ ^^e associated pair 
basis the diagonal matrix elements of the Hamiltonian 
1-L correspond to the 'bare' binding energies of the pair 
states; i.e. the binding energy in the abscence of coupling 
between the channels by V{r). Restricting ourselves, for 
purposes of introduction, to a single open channel and to 



the least bound states in the interaction potentials, Hpp 
is a single matrix element on the diagonal of H, corre- 
sponding to the hare binding energy of the least bound 
state of the entrance channel, ep = —h^H?p/2ii. The en- 
ergy ep can be readily calculated by projecting the pair 
state on the spin basis {|cr) = \SMsiJial~^^)}i and is given 
by 

^P=^4^ Yl {SMsl^al^(3\famf^Jf3mf^f (18) 

In the following sections we will show that the width AB 
of the resonance is a function of the bare binding energy 
ep of the entrance channel and a single matrix element 
of ^PQ^ denoted by /C, representing the coupling of the 
level ep to the resonant level in Q space. We will show 
that the width is given by the expression 



jlrel^B 



K^ 



2ai,gnp\ep\' 



(19) 



Hence, once the ABM parameters are known, the width 
of the resonances follows with an additional unitary 
transformation of the ABM matrix to obtain the cou- 
pling coefficient /C. In view of the orthogonality of the 
subspaces V and Q, the submatrix Hqq^ corresponding 
to all closed channels, can be diagonalized, leaving Hpp 
unaffected but changing the Hpq and Hqp submatri- 
ces. In diagonalized form the Hqq submatrix contains 
the energies eg of all bound levels in Q space and in- 
cludes the coupling of all channels except the coupling 
to V space. This transformation allows to identify the 
resonant bound state and the corresponding off-diagonal 
matrix element /C in Hpg , which is a measure for the res- 
onance width. Thus we can obtain the coupling between 
the open and closed channels without the actual use of 
continuum states. 

In Section |IVB| we present the Feshbach theory tai- 
lored to suit the ABM. We give a detailed description 
of the resonant coupling, and demonstrate with a two- 
channel model how the ABM bound state energy E^j com- 
pares to the associated 7^-space bare energy ep, and to 
the dressed level in the entrance channel from which one 

we gen- 



can deduce the resonance width. In Section lIVC 



eralize the results such that the width of the Feshbach 
resonances can be obtained for the general multi-channel 
case. For a more thorough treatment of the Feshbach 
formalism we refer the reader to [25l [26] and for its ap- 
plication to cold atom scattering e.g. [TO] , 



B. Tailored Feshbach theory 

By introducing the projector operators P and Q, which 
project onto the subspaces V and Q, respectively, the 
two-body Schrodinger equation can be split into a set of 
coupled equations [10 



{E-npp)\^p) = HpqI^q), 



(20) 



iE-nQQ)\<bQ) = HqpI^p), 



(21) is reduced to only one term. Therefore, the complex en- 



where \^p) = P|^), |^g) = Q|^), Hpf = PUP, 
HpQ = PHQ^ etc. Within Q space the Hamiltonian 
Hqq is diagonal with eigenstates \(J)q) corresponding to 
the two-body bound state with eigenvalues eg. The en- 
ergy E = h'^k'^/2/jj is defined with respect to the open 
channel dissociation threshold. 

We consider one open channel and assume that near 
a resonance it couples to a single closed channel. This 
allows us to write the 5* matrix of the effective problem 
in V space as [10] 



+ \|2 



.,.,^«., :-...KM^ 



(22) 



where |^p) are scattering eigenstates of Hpp, Sp{k) is 
the direct scattering matrix describing the scattering pro- 
cess in V space in the absence of coupling to Q space. 

The complex energy shift A{E) describes the dressing 
of the bare bound state \(J)q) by the coupling to the V 
space and is represented by 



A{E) = {<Pq\Hqp 



1 



E+ -U 



pp 



-npQ\<PQ), (23) 



where E~^ = E -\- iS with 5 approaching zero from 
positive values. Usually the open channel propagator 
[E~^ — 1-Lpp]~^ is expanded to a complete set of eigen- 
states of Hpp^ where the dominant contribution comes 
from scattering states. To circumvent the use of scatter- 
ing states we expand the propagator to Gamow resonance 
states, which leads to a Mittag-Leffler expansion [27j 



E+ - Upp 






E 



\^n){^^\ 



^ kn\k kn) 



(24) 



where n runs over all poles of the Sp matrix. The 
Gamow state \Qri) is an eigenstate of Hpp with eigen- 
value ep^ = h'^k'^/{2fi). Correspondingly, the dual state 
|n^) = |r^n)*7 is an eigenstate of Wpp with eigenvalue 
(ep^)*. Using these dual states, the Gamow states form a 
biorthogonal set such that {rt^\rtn') = ^nn' • For bound- 
state poles kn = ii^m where Kn > ^^ Gamow states cor- 
respond to properly normalized bound states. 

We assume the scattering in the open channel is dom- 
inated by a single bound state {kn = inp). This allows 
us to write the direct scattering matrix in Eq. (22) as 



Sp{k) 



_ g-2i/cabg 



-2ikaL f^P 



ik 



tvp + ik 



(25) 



where abg is the open channel scattering length, and the 
P-channel background scattering length a^ is on the or- 
der of the range of the interaction potential a^^ ^ tq. 
Since we only have to consider one bound state pole (with 
energy ep) in P space, the Mittag-Leffler series Eq. (|24|) 



ergy shift Eq. (23) reduces to 



A{E) 



M 



-zA 



h"^ Kp{k — IKp) 



(26) 



where A = ((/)q|Hqp|1^p)(1^p|Hpq|(/)q) is a positive 
constant. The coupling matrix element between open- 
channel bound state and the closed-channel bound state 
responsible for the Feshbach resonance is related to A. 

The complex energy shift can be decomposed into a 
real and imaginary part such that A{E) = A^esiE) — 
^r{E). For energies E > the unperturbed bound state 
becomes a quasi-bound state: its energy undergoes a shift 
Ares and acquires a finite width F. For energies below the 
open-channel threshold, i.e. E < 0^ A{E) is purely real 
and r{E) = 0. In the low-energy limit. A: ^ 0, Eq. ([26| 
reduces to 



A{E) = A - iCk, 



(27) 



where C is a constant characterizing the coupling 
strength between V and Q space [10 , given by C = 
A(2A^p|ep|)-^ and A = A(2|ep|)-^ Note that if the 
direct interaction is resonant, |abg| ^ ro, the energy de- 
pendence of the complex energy shift is given by [28] 
A{E) = A — iCk{l + ika^)~^ where a^ = K,p^^ yielding 
an energy dependence of the resonance shift. 
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Figure 6: (Color online) Illustration of the threshold behavior 
in a (fictions) two-channel version of the dressed ABM. The 
threshold behavior is determined by the coupling between the 
least bound level in the open channel in V space with the reso- 
nant bound level in Q space. The uncoupled levels are shown 
as the blue (ep) and red (eg) dash-dotted lines, with eg cross- 
ing the threshold at Bo- The solid black lines represent the 
dressed levels, with the upper branch crossing the threshold 
at B{). Near the threshold, the dressed level shows the char- 
acteristic quadratic dependence on (B — Bq) (see inset). For 
pure ABM levels (dotted gray) no threshold effects occur and 
the coupled bound state crosses the threshold at Bq. 

Since we consider one open channel, the (elastic) 5*- 
matrix element can be written as e^^^^^\ where 5{k) is 



10 



the scattering phase shift. The scattering length, defined 
as the limit a = — tan5(/c)//c, {k -^ 0), is found to be 
Eq. ([3| which shows the well known dispersive behavior. 
The direct scattering process is described by the scatter- 
ing length ttbg = a^ -\- a^ . At magnetic field value ^o, 
where the dressed bound state crosses the threshold of 
the entrance channel, the scattering length has a singu- 
larity. 

The dressed state can be considered as a (quasi-) 
bound state of the total scattering system. The energy of 
these states is obtained by finding the poles of the total 
S matrix Eq. p2|). This results in solving 



(28) 



{k - inp) {E-CQ- A{E)) = 0, 



for k. Due to the underlying assumptions, this equation 
is only valid for energies around threshold where the open 
and closed channel poles dominate. Near threshold, the 
shifted energy of the uncoupled molecular state, eg + A, 
can be approximated by iirei^B — Bq). This allows to 
solve Eq. (28) for E and we readily obtain 



E 



2\ep\y^firei{B-Bo) 



(29) 



retrieving the characteristic quadratic threshold behavior 
of the dressed level as a function of {B — Bq). The energy 
dependence of molecular state close to resonance is also 
given hy E = —h^/{2iio?) this allows us to express the 
field width of the resonance as A5 = C{a\y^iirei)~^ • 

We apply the above Feshbach theory to a (fictitious) 
two-channel version of ABM, and the results are shown in 
Fig. [6J This two-channel system is represented by a sym- 
metric 2x2 Hamiltonian matrix, where there is only one 
open and one closed channel. The open and closed chan- 
nel binding energies ep resp. eg are given by the diago- 
nal matrix elements, while the coupling is represented by 
the (identical) off-diagonal matrix elements. The closed 
channel bound state is made linearly dependent on the 
magnetic field, while the coupling is taken constant. In 
addition to ep and eg, we plot the corresponding ABM 
solution, which in this case is equivalent to a typical two- 
level avoided crossing solution. The figure now nicely 
illustrates the evolution from ABM to the dressed ABM 
approach, where the latter solutions are found from the 
two physical solutions of Eq. (28), which are also plotted. 



Since the dressed ABM solutions account for threshold 
effects, they show the characteristic quadratic bending 
towards threshold as a function of magnetic field. From 
this curvature the resonance width can be deduced. 



C. The dressed Asymptotic Bound state Model 

To illustrate the presented model for a realistic case, 
we will discuss the ^Li-^^K system prepared in the 
|/Limj^,,/Kmj^) = |l/2,+l/2,9/2,-7/2) two-body hy- 
perfine state as an example throughout this section. This 
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Figure 7: (Color online) Dressed bound states for ^Li-^'^K for 
Mf — —3 (black lines, see also Table [l|. The uncoupled Q 
and V bound states (I-Lpq = 'Hqp = 0) are represented by 
the dot-dashed lines (red and blue respectively). The gray 
shaded area is shown in detail in Fig. |8] 



particular mixture is the energetically lowest spin com- 
bination of the Mp = —3 manifold, allowing to consider 
only one open channel. We note that the model can be 
utilized to cases containing more open channels. 

In order to calculate the width of a Feshbach resonance 



using the method presented in Sect. |IVB| three quanti- 
ties are required: the binding energy of the open channel 
ep, of the closed channel responsible for the Feshbach res- 
onance eg, and the coupling matrix element /C between 
the two channels. In the following we will describe how 
to obtain these quantities from the ABM by two simple 
basis transformations. 

For ultracold collisions the hyperfine and Zeeman in- 
teractions determine the threshold of the various chan- 
nels and thus the partitioning of the Hilbert space into 
subspaces V and Q, and therefore a natural basis for our 
tailored Feshbach formalism consists of the eigenstates of 
W^^. Experimentally a system is prepared in an eigen- 
state \al3) of the internal Hamiltonian W^^, which we 



refer to as the entrance channel (cf. Section II A). Per- 
forming a basis transformation from the \a) basis to the 
pair basis \aP) allows us to identify the open and closed 
channel subspace. The open channel has the same spin- 
structure as the entrance channel. 

We now perform a second basis transformation which 
diagonalizes within Q space without affecting V space. 
We obtain the eigenstates of Hqq and are able to iden- 
tify the bound state responsible for a particular Fesh- 
bach resonance. The bare bound states of Q space 
are defined as {|^gi), |</>g2)7 • • •} with binding energies 
{eg^, eg2, . . .}. For the one dimensional V space, which 
is unaltered by this transformation, the bare bound state 
\^p) of Hpp is readily identified with binding energy 
ep. In the basis of eigenstates of Hpp and Hqq we 
easily find the coupling matrix elements between the i- 
th Q space bound state and the open channel bound 
state ((/>gjHgp|r^p). This gives the coupling constant 
Ai = {<pQ,\nQp\Qp){Q^\npQ\<PQ,) = IC that deter- 
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mines the resonance field Bq by solving Eq. ( 28 ) at thres- 
hold, 



^2 



eg^ep 



(30) 



The field width of this Feshbach resonance is proportional 
to the magnetic field difference between the crossings of 
the dressed (Bq) and uncoupled Q bound states (Bq) 
with threshold since 
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Figure 8: (Color online) A zoom of the dressed ABM as shown 
in figure|7| The dressed molecular states are shown near thres- 
hold (black) . The field width of a resonance is related to the 
magnetic field difference of where the dressed and uncoupled 
Q bound state cross the threshold. 



Table I: The positions of all experimentally observed s- 
wave Feshbach resonances of ^Li— ^°K. Column 2 gives the 
^Li (mfj^.) and ^"^K (rnfj^) hyperfine states. For all reso- 
nances fLi = 1/2 and fx = 9/2. Note that the experimental 
width of the loss feature ABexp is not the same as the field 
width A^ of the scattering length singularity. Feshbach res- 
onance positions Bq and widths AB for ^Li-^°K as obtained 
by the dressed ABM, obtained by minimizing x^ • The last 
two columns show the results of full coupled channels (CC) 
calculations. All magnetic fields are given in Gauss. The ex- 
perimental and CC values for Mf < and Mf > are taken 
from Ref. 8 and [9] respectively. The resonances marked 
with * have also been studied in Refs. [291130] . 
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e^'^. Therefore, we have again preformed a y^ analysis 
and we found new values of the binding energies e^ /h = 
713 MHz and e^/h = 425 MHz, which yields a lower x^ 
minimum as compared to the regular ABM calculation. 
The obtained value of AB generally underestimates 
the field width of a resonance. This originates from the 
fact that only the dominant bound state pole correspond- 
ing to a^ has been taken into account. By including the 
pole of the dominant virtual state in the Mittag-Leffier 
expansion, the coupling between the open and closed 
channel will increase, hence, AB will increase. 



V. SUMMARY AND CONCLUSION 



We illustrate the dressed ABM for ^Li-^^^K in Figs, (t] 
and |8J for Mp = — 3. To demonstrate the effect of 
^PQi we plotted for comparison both the uncoupled and 
dressed bound states [35 . Details of near-threshold be- 
havior (gray shaded area in Fig.[7| are shown in Fig.jsjto- 
gether with the obtained scattering length. We solved the 
pole equation of the total S'-matrix Eq. ( 28 ) for each Q- 



state and plotted only the physical solutions which cause 
Feshbach resonances. The dressed bound states show the 
characteristic quadratic bending near the threshold. We 
have used Cq to determine tq (~ a^ ) from Eq. 16 

Table U summarizes the results of the dressea ABM 
for the ^Li-^^K mixture. Note that the position of the 
Feshbach resonances will be slightly different compared 
to the results from the regular ABM, for equal values of 



We have presented a model to accurately describe 
Feshbach resonances. The model allows for fast and ac- 
curate prediction of resonance positions and widths with 
very little experimental input. The reduction of the ba- 
sis to only a few states allows to describe Feshbach res- 
onances in a large variety of systems without accurate 
knowledge of scattering potentials. 

Using the ABM in combination with the accumulated 
phase method allows to describe Feshbach resonances in 
alkali systems with a large degree of accuracy, using only 
three input parameters. Additionally, the fast compu- 
tational time of the model allows to map all available 
Feshbach resonances in a system and select the optimal 
resonance required to perform a certain experiment. For 
the ^Li-^^K mixture we have utilized this ability to find 
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a broad resonance as presented in Ref. [9]. In addition, 
locating e.g. overlapping resonances in multi-component 
(spin) mixtures can be performed easily well using the 
ABM. 

An additional important feature is that the model can 
be stepwise extended to include more phenomena allow- 
ing to describe more complex systems. For example, a 
possible extension would be to include the contribution 
of the dominant virtual state in the Mittag-Leffler expan- 
sion, this would allow for an accurate description of the 
resonance widths for systems with a large and negative 
ttbg. Additionally, including the dipole-dipole interaction 
allows to describe systems where Feshbach resonances oc- 
cur due to dipole-dipole coupling. Finally, it has already 
been shown by Tsherbul, et al. [ST that the ABM can 
be succesfully extended by including coupling to bound 



states by means of an externally applied radio-frequency 
field. 

The approach as described in the ABM is in principle 
not limited to two-body systems. Magnetic field induced 
resonances in e.g. dimer-dimer scattering have already 
been experimentally observed [32 . For few-body systems 
an approach without having to solve the coupled radial 
Schrodinger equations is very favorable. This large vari- 
ety of unexplored features illustrate the richness of the 
model. 

This work is part of the research program on Quan- 
tum Gases of the Stichting voor Fundamenteel Onder- 
zoek der Materie (FOM), which is financially supported 
by the Nederlandse Organisatie voor Wetenschappelijk 
Onderzoek (NWO). 
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